# import the fdc dataset --------------------------------------------------

df_fdc <- read.csv("datasets/data_fdc.csv")

# import the map of turkey ------------------------------------------------

tr_map <- readRDS("datasets/data_map.RDS")

# plot the figure ---------------------------------------------------------

figure_2 <- tr_map +
  geom_point(data = subset(df_fdc, casualties > 0),
             aes(x = longitude, y = latitude, size = casualties),
             color = "red", shape = 1, stroke = 1.2) +
  coord_sf() +
  labs(x = NULL, y = NULL) +
  theme_minimal() +
  theme(panel.grid = element_blank(),
        axis.text = element_blank(),
        legend.position = "bottom") +
  scale_size_continuous(name = "Casualties", range = c(2, 6))

# display the figure ------------------------------------------------------

print(figure_2)

# save the figure ---------------------------------------------------------

ggsave(plot = figure_2, filename = "figures/figure_2.pdf", dpi = 1000,
       width = 9, height = 4, units = "in")